Statistics of resonances in a one-dimensional chain: 
a weak disorder limit 



Vinayak^'^ 

^Institute de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, 
Cuernavaca, Mexico 

^Department of Physics, Technion-Israel, Institute of Technology, Haifa 32000, Israel 

Abstract. We study statistics of resonances in a one-dimensional disordered chain 
coupled to an outer world simulated by a perfect lead. We consider a limiting case 
for weak disorder and derive some results which are new in these studies. The main 
focus of the present study is to describe statistics of the scattered complex energies. 
We derive compact analytic statistical results for long chains. A comparison of these 
results has been found to be in good agreement with numerical simulations. 



PACS numbers: 03.65.Yz, 03.65.Nk, 72.15.Rn 



Statistics of resonances in a one- dimensional chain: a weak disorder limit 



2 



I. Introduction 

Resonant phenomena have received much attention in atomic and nuclear physics and 
more recently in chaotic and disordered systems [H O [3l HI O [6l [7] . Complex energies, 
Ea = Ea — ^Ta, which correspond to poles of the scattering matrix on the unphysical 
sheet, characterize resonances [S]. Resonances correspond to the long-lived quasi- 
stationary states which eventually decay to continuum while distribution of resonance 
widths, -P(r), determines decay of the corresponding survival probability with time. 

In recent years, P(T) has been a subject of investigations jH [2] for a, simple but 
much studied, discrete tight-binding one dimensional random chain which is coupled 
to a perfect lead at one side. A numerical study [2] shows that in a broad range 
of r, P(r) ~ r~'^, where the exponent 7 is very close to 1. Intuitively the l/F 
behaviour can be deduced by assuming a uniform distribution for the localization centers 
of exponentially localized states |3]. However, from analytic point of view one usually 
considers an infinitely long chain in which case the average density of resonances (DOR) 
has a well defined limit. For a finite size system, the difference between the DOR 
and -P(r) is the normalization by the system size [H |5]. Recently, Kunz and Shapiro 
have derived analytic expression of the DOR for a semi- infinite disordered chain [5]. 
They have obtained an exact integral representation of the DOR which is valid for 
arbitrary lead-chain coupling strength. This has been further simplified for small lead- 
chain coupling strength where a universal scaling formula is found. In this limit they 
have proved the l/F-behavior of the DOR [U E]. Besides, for the continuous limit of 
this model an integral representation of DOR has been obtained |6]. 

Kunz and Shapiro's work has established a universal l/F law for arbitrary strength 
of disorder in a semi-infinite chain. Numerically one can verify l/F law of the DOR, 
similar to what has been done by Terraneo and Guarnery [2] in finite samples for -P(F). 
Such verifications require the localization lengths to be much smaller than the size of 
the sample. In case of weak disorder an analytic result for the localization lengths is 
particularly useful. It comes from a second order perturbation theory. It states that the 
localization length is maximum near the middle of the energy band and is proportional 
to W'"^ where W is the width of the disorder [9l [TOl [11] . On the other hand, this result 
also leads to an interesting limiting situation where the localization lengths are much 
longer than the sample size. This is what we refer to as a weak disorder limit in this 
paper. This limit has scarcely been studied hitherto although it is relevant in the study 
of localization through resonances. Besides, there has been a believe for some sort of 
universality in the weak disorder limit. In this paper we address to this limit and derive 
analytic results which describe the statistics of resonances. Our work probes a fresh 
area and studies a weak disorder limit which has never been addressed before. 

For open systems, instead of studying the scattering matrix in a complex plane we 
follow an alternative approach where one solves the Schrodinger equation by describing a 
particle ejected from the system or equivalently with a boundary condition of outgoing 
waves (Siegert boundary condition [I2]). In this approach one naturally turns up to 
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Figure 1. A one-dimensional disordered chain with A'' sites, represented in the figure 
by black dots, is coupled to a lead. Open circles represent sites of the lead. The 
outgoing plane wave is shown by the arrow where < 5R{fc} < tt and < 0, so 

that it propagates left in the lead and its amplitude grows in the lead. 



a problem of solving a non-Hermitian effective Hamiltonian which admits complex 
eigenvalues |21 SI El [13] • For details of such non-Hermitian effective Hamiltonians, 
we refer to a recent study [H] and references therein. 

We derive the statistics which describe scattered complex energies of disordered 
chain around those regular ones which correspond to an open chain without any disorder 
(clean chain). For instance, we derive average of square of the absolute values of the 
shifts in complex energies from the regular ones over all realizations of the set of random 
site energies. Similarly we obtain results for the statistics of real and imaginary parts 
of those shifts. These results lead to compact expressions for long chains. To show 
the generality of our approach we also derive these results for the so-called parametric 
resonances which have been particularly useful in numerical studies [2] . Finally, we give 
numerical verifications of our analytic results. 

The paper is organized as follows. Although the system and its effective 
Hamiltonian have been nicely explained earlier in [21 HI |5], for the sake of completeness 
of this paper we will describe these briefly in section H. In the same section we will also 
describe the exact and the parametric resonances. In Sec. HI we will derive result for 
resonances in an open-clean chain of finite length, in terms of a polynomial equation. 
For long chains, we will solve this polynomial equation in the leading order of the inverse 
of the length. In Sec. IV we will use the perturbation theory to obtain the first and 
the second order corrections in the complex energies for a weak disorder. In Sec. V 
we will calculate statistics of the scattered complex energies. In the same section we 
will simplify our results for long chains and obtain compact expressions. In Sec. VI we 
will briefly discuss about the numerical methods to calculate complex energies of non- 
Hermitian effective Hamiltonians and numerically verify our analytical results. This will 
be followed by the conclusion in Sec. VII. 

II. Model and Its Effective Hamiltonian 

A discrete tight-binding one dimensional chain of length N (shown by positive integers, 
n = 1, 2, A^, used for indexing the sites of the chain in Fig. [1]) is connected to an 
outer world (represented by a perfect lead which sites are shown by a zero and negative 
integers, n = 0, —1, —2, ...). Each site of the chain has the site energy e„ where e„ are 
statistically independent random variables chosen from some symmetric distribution. 
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Each nearest neighbor site of the chain as well as of the lead is coupled by a hopping 
amplitude t. The hopping amplitude for the pair n = and n = 1 is t' which takes 
values from t' = (closed chain) to t' = t (fully coupled chain). With this hopping, 
a particle, which is initially located somewhere in the chain, eventually escapes to the 
outer world. 

Now we write down the Schrodinger equation for the entire system. 



- ti)n+l - tiJn-l = S'ipn, for H < 0, (l) 

- t'l/ji = ^Vo, for n = 0, (2) 

-t'ipo - til)2 + ^ii>i = Si^i, for n = 1, (3) 

-t^„_l - tV„+i + eni^n = ^i^n, foY 2 < U < N . (4) 



In order to avoid cluttering of notations we always represent quantities corresponding to 
disordered system by script letters while quantities for the clean system are represented 
in usual math notations. Tilde is used to discriminate the open system case from the 
closed one. Equation ([T]) is for the lead where e„ = 0. Equations (O [3]) describe the 
lead-chain coupling and Eq. (jlj) is for the chain. As in, [5] we solve Eqs. ([iMj) with 
a boundary condition of an outgoing plane wave in the lead, i.e., ipn<o oc exp{—ikn) 
where < ^{k} < tc and '^{k} < 0. The condition on ^{k} ensures that the outgoing 
wave propagates to left in the lead. The condition on '^{k} is considered so that the 
amplitude of the resonance wave function grows in the lead. It comes from Eq. ([T]) that 
the complex energy S is related to the complex wave vector k via the dispersion relation 
S = — 2tcos(fc). Now we eliminate all ipn for ^ < 1 from Eqs. ([HH]) and obtain 

-t^Jn+l - t^Jn-l + in^n = S^n, (5) 

where 

en = en-tT] exp{ik)6ni, (6) 

for n = 1, 2, A^. The parameter rj = {f /tY measures the coupling strength to the 
outside world. 

An effective Hamiltonian defined by the Eq. ([5]) is non- Hermit ian. For instance, if 
T-L is the N x N tridiagonal Hermitian matrix which represents the Hamiltonian of the 
closed-disordered chain then one may write the effective Hamiltonian, 7/, as 

n = n-tr]XCk)P. (7) 

Here P = |1)(1| is the projection for site n = 1 and A = exp{ik). The above non- 
Hermitian effective Hamiltonian has been first obtained by Terraneo and Guarnery [2]. 
The underlying result here is that the same relation (JTj) is valid for any Hermitian 
H representing a (closed) quantum system [H] which has A^- dimensional state space. 
Resonances are characterized by the complex eigenvalues, Sa, of Ti. 

Note here dependency of "H on the complex wave vector k which is related to the 
complex energies via the dispersion relation mentioned above - this is not a standard 
eigenvalue problem. To standardize this problem "parametric resonances" are often 
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used as an alternative. In this approach the dependence of A on A; is typically neglected, 
reducing thereby the problem of finding the eigenvalues of the effective Hamiltonian 
at chosen value of k. As expected, parametric resonances yield approximate statistical 
results which are close to those for the exact resonances in strongly localized regime [2]. 
Parametric resonances depend on a chosen parameter, for instance let k = ko and we 
fix it in the middle of the energy band, k^ = tt/2. Writing explicitly 

\exT)(ik), for exact resonances, , ^ 

m = . ' . (8) 

tor parametric resonances. 

From now on we set the energy scale by taking t = 1, denoting the complex variable 
S/thyZ. We denote the Hamiltonian matrix representing the closed-clean chain by H. 
It differs from "H only at the diagonal as, for the clean chain case, all the site energies 
are zero. Calculation of the eigenvalues of if is a standard exercise where one derives 
Za = -2 cos[a7r/(A^ + 1)] for a = 1, A^. 

Before going into a detail treatment to the problem, we should first sketch the 
outline of our approach. We are interested in a weak disorder regime. Since our approach 
rely on perturbation theory, we need complex energies of open-clean chain, i.e., the Zq,s. 
So we will begin with calculating the resonances for open-clean chain of finite length. 
Then we will do the perturbation series expansion up to the second order of strength of 
the disorder. This will be followed by the derivation of the statistical results. Finally, 
we will consider the large- iV limit of these results. 



III. Open-clean Chain 

We begin with defining the resolvent Q{z) = {z — l-L)^^ . Using Eq. ([7]) we may also 
write 

g{z) = {z-H + r]XP)-\ (9) 
For the open-clean chain we define the resolvent 
G{z) = {z-H + r]XP)-^ 

= {1 + 7]XGP)-^G, (10) 

where we have introduced G{z) = {z — H)^^ as the resolvent for the "unperturbed" 
closed-clean chain. Resonances correspond to the singularities of the matrix Grnn{z), or 
to the roots of the secular equation 

F{z) = = l + r]XG^,iz), (11) 

where Gu is the {1, 1} element of the matrix G in site representation. (G„m(-2) = 
{n\{z- H)~^\m).) 

To obtain Gu for finite A^, we use the ordinary difference equation (ODE), 



(12) 



Statistics of resonances in a one- dimensional chain: a weak disorder limit 



6 




Figure 2. Comparison of the result ((26|) (pluses) with the numerical solution of the 
polynomial equation (I20p (circles, squares and diamonds) for the exact resonances 
where 77 = 0.5, 0.81 and 0.99. We have considered N — 100. 



with the boundary conditions ipo = ipN+i = 0. This equation is obtained from Eq. 
([5]) by setting all e„ = 0. Next we consider Un{z) and Vn{z) to be the two linearly 
independent functions which satisfy the ODE 

M„+l + + ZM„ = 0, (13) 

Vn+l+ Vn-l + ZVn =0, (14) 

where uq = vn+i = 0. Since norm of Un, Vn is arbitrary, we fix -ui = wat = 1. Further 
we claim that the resolvent is given by 

_ UnVmQ{rn - n) + UmVnOin - m) 

Here Q{n) is the unit-step function and Wn = UnVn-i — Un-iVn is the Wronskian. Using 
Eqs. ( fT3l [T4l) it is straight forward to see that the Wronskian is independent of n. One 
can also check that 

Gn+lm ~\~ Gn—lm ~l~ ^ Gnm ^nm- (-^6) 

We now set u„ = vjy+i^n to match the initial value problem f lT^ [T^ to the boundary 
value problem (fT2|) . We find 

Gn = (17) 
un+i 

The ODE f lT3|) is satisfied by the Chebyshev polynomial of the second kind, Um{—z/2), 
defined as 

, . sinffm + 1) cos~"'^(a;)l 
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Figure 3. Repeated on the same pattern of Fig. [2] but for parametric resonances. 

for Uq{x) = 1 and Ui{x) = 2x. Since we have fixed Ui = 1, therefore = f/„_i, thus 
we can write Eq. ( !T7l) as 

Un^i{-z/2) sm[Nk] 



Gu 



(19) 



Un{-z/2) sm[{N + l)ky 

Here the last equahty follows from the energy dispersion relation. Using Eq. ( !T9|) in 
Eq. ( ITTl) . we end up with an algebraic equation 

UN-ii-z/2) 



F{z) = = 1 - r/A- 



(20) 



Un{-zI2) ■ 

Zeros of F{z) are the roots of a polynomial of order A^. For exact resonances Eq. 
( l20l) can be easily transformed into 

(2iv+i) _ 3{z)-^ - ri 3{z) 



[a(^)]^ 



1 — 1] 



(21) 



where 

a(z) = -exp[iA;(z)]. (22) 

In order to solve Eq. (12T]) . we propose an ansatz assuming that opening of the system 
at one end causes 0{N~^) complex corrections to the /Cq,'s. Let 



iV 



(23) 



where is a complex quantity and ka = air /[N + 1). Inserting this ansatz into Eqs. 
(EH [22D we obtain 



i ( 1 

= - — In [n(/c«; r/)] + O ( — 



(24) 
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Figure 4. Scatter plot for exact resonances where N — 100, rj = 0.81. Dense points 
in the graph represent exact resonances in the disordered chain for 2500 reahzations 
where W = 0.015. These are scattered around dots which represent exact resonances 
in the clean chain. 



where 

Now, up to 0{N~^), Za may be written as 

Za = —2 cos(A;q,) 



?7 e 



2^ k(y 



i sin(fco 



ln(^]). 



(25) 



(26) 



The same result can be obtained for the parametric resonances, after repeating the 
similar steps, but with different Q: 



irj e 



1 — trje 

One should bear in mind that there is no resonance for rj 



(27) 



1, as the system is fully 
coupled to the lead. However, for parametric resonances, one artificially gets resonances 
even when rj = 1. Note that the result fl26l) is symmetric about the imaginary axis for 
both cases. In Fig. [2] and Fig. [3] we compare the numerical solutions of the polynomial 
equation ([20]) with our results ([261 ^ IIZD, for = 100, r] = 0.50, 0.81 and 0.99 and 
= 100, respectively for exact and parametric resonances. Eq. (1201) has been solved 
by using the Newton's method with the initial guess = k^. These figures show that 
our result fl26|) is close to the numerical solution. The agreement gets better as — 1 
(not shown here separately). However, the ansatz ( 123|) is not valid near the band edges. 
Moreover, the agreement fails for parametric resonances near the middle of the band as 
?7 — )■ 1; see Fig. [3] for 77 = 0.99. 



Statistics of resonances in a one- dimensional chain: a weak disorder limit 



9 




S -0.015 



-0.025 



Re{Z } 

Figure 5. Scatter plot for parametric resonances wliere N — 100, rj = 0.81 and 
W = 0.015, for 5000 realizations. As in FigHl here also dots represent the clean chain 
and points represent the disordered chain. 



IV. The Weak disorder limit 



In the next stage of the problem we switch on a very weak disorder in the chain. From 
a second order perturbation theory we know that for a disordered infinitely long chain 
the localization length, ^(-E), is maximum at the middle of the band. For small W it is 
given by [101 E] 

24(4t2 - E^) 



(28) 



However, the exact result shows a small deviation at 



implying thereby, ^(0) 

the band center due to the breakdown of the second-order perturbation theory [11]. We 
consider a limiting situation when ^(0)/A^ >> 1. For instance, in Fig. [Hand in Fig. 
El we show the scatter plot (JR.{Zo,} vs '^{Za}) for exact and parametric resonances 
respectively. In both cases we have considered N = 100, rj = 0.81 and W = 0.015 so 
that ^(0) >> N. As seen in these figures, complex energies of the disordered chain are 
scattered around the ZaS. 

We now calculate the corrections to Za for such weak disorder case. It is suggestive 
here to deal with the self-energy. Let iSi(e2, eN', z) be the self-energy for the first site, 
defined via 

^ii(^) = . V (29) 

Here {e} denotes the set e2,...,eAr and Qu is the {1, 1} element of the resolvent 
Q{z) = {z — "H)"^, defined for the Hermitian matrix "H. For the later convenience 
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we write 



n = H + w, 



(30) 



where W = J2'e=i ^t^t ^^'^ Pi= K ) ( ^ I is the projection for the £'th site. 

In the rest of the paper we will work out results only for the exact resonances. For 
the parametric resonance theses results can be carried out following similar steps, so we 
skip all the intermediate steps merely by stating the result at the end. 

As before in Eq. ( ITTl) . for disordered chain, resonances correspond to the roots of 
the secular equation 

jr(^z)=f) = z-ei-Si{{t}]z) + \r]. (31) 

Preserving Za as the roots of Eq. ( fTTl) . we define as the roots of Eq. ( 13T]) . Now 
we expand the roots = 5^ + {5iZa) + (52-2^), assuming that {5iZa) are linear while 
(52^0,) are quadratic in the , for j = 1,...,A^. Then for Si{{e}] Za), up to C({e}^), 
we get 

N 



n=2 V 



-0,Z=Za 



+ (SlZa + ^2^0 



N 



n,m=2 



-(5l4 



dSr{{e};z) 
dz 

d'S,i{e};z] 

d^S,{{e}-z) 
d^z 



{.}=o„ 



{e}=0,Z=2a 



{e} = 0,Z = Za 



We will use this expansion in Eq. ( |3T]) . Before that we evaluate 
'dS^{{t}-z)\ _ d 1 

/ {e}=0,z=Zc 



dz 



and. 



dS,i{ey,z) 



der 



dz Gii{z) 



(32) 



(33) 



{e}=0,Z = Za 



{Guy der. 



{e} = 0,Z = Za 



(34) 

for n > 2. These equalities come from Eq. ( !29l) . Finally, we calculate derivatives of 
Qii, at {e} = and z = Za with respect to {e} by using Eq. ( l30i) for the Born-series 
expansion of Q{z). We find 



der. 



— GlnGnl 

{e}=0,Z = Za 

Grouping all these, for the first order corrections, we obtain 



(35) 



{SiZa] 



ei + En= 


GlnGnl 

[Giif 


Z — Zen 


d 1 


^ iriexp[ika] 
z=z^ 2sm{ka) 


dz Guiz) 



(36) 
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Similarly for the second order corrections we get 



X 



N 

E 

n,m=2 



d 



GlnGnmGml [GlnG 



Im I 



[Gil? 



d^z Gu 



dz Gii{z) 



+ 



if] ey.\;)[ika 
2 sin(A;Q,) 



[Giif 

d"^ exp{ik{z) 
d^ 

-1 



{e}=0,Z = Za 



(37) 



Note that (SiZa) and (52^'^) have been obtained in terms of the resolvent of the 
closed-clean chain which we already know in terms of Chebyshev polynomials; see Eq. 
( !T5|l and the relation between u„ and f„ with Chebyshev polynomials. 



V. Statistics of the Scattered Complex Energies 

We are interested in the statistics of the scattered complex energies. For instance, using 
the first order result (1361) of the perturbation theory, we calculate average of square of 
absolute shift in complex energies defined as, (|(AZq,)P) = (|(^q, — 5a)P)- The angular 
brackets are used here to represent the averaging over many realizations of set of all 
random site energies {e„}. This quantity gives a statistical account for the scattered 
complex energies. We also calculate {(JR,{AZa}Y) and ( (53{AZa})^ ), viz, average of 
square of the real and the imaginary part of the shift {Za — Za), respectively. To obtain 
the latter quantities we need first to calculate ( (AZa) ^) and ( [{AZa)*? ), since 



i^{AZ^}) 



[Q{AZ^}Y 



2_ (AZ,)2 + [(AZ,)*]2 + 2(|AZ„)|2) 



4 



(AZ,)2 + [(AZ«)*]2-2(|AZ«)P) 



(38) 



(39) 



Here we have used {*} to represent the complex conjugate (c.c). 

For all these three statistics we simplify {6iZa), given in Eq. ( I36ll . in terms of 
Chebyshev polynomials as 

N 



(SiZ^) 



£-4 



n=l 



X 



?7jv-i(5«/2)T^+i(5j2)-iV 



ir]exp[ika] sm{ko,)[UN-i{za/2)f 



-1 



(40) 
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where Tm{z) = cos[m cos~^(2)] is the Chebyshev polynomial of the first kind. Further 
simplifications occur when these polynomials are expressed in their trigonometric forms. 
For instance, let's calculate |(AZq,)P, with 5q,/2 = cos(6'q,) where 9a = n — ka- We obtain 

4 \D{zaW ■ ^ ' 

Here n' and m' are respectively N + 1 — n and N + 1 — m, and D{za) is simply the 
quantity in the second bracket of Eq. (HOl) . Averaging releases one of the summation 
as the e^s are statistically independent-identically-distributed (i.i.d.) random variables. 
We simply have 

(l(A^.)r) = a^Eiil^!!^^^^^ (42) 

where cr^ is variance of the e^s. 

Summation in the above equality can be performed by using trigonometric 
identities. For instance, we first write 

4sin2(n^) sm^{n9*) = 1 - cos(2n0) - cos{2n9*) 

^ cos(4 n ^{9}) + cos(4 z n S{g}) ^^^^ 



and we use the summation formula 

N 

I einll /V -i- I /y\Hl 

(44) 



Yl cos(ri^) = i 

n=l 



sin[(iV+l/2)^^] 



sm{9/2) 



It turns out after some trigonometry that one can write the summation in a closed form. 
We find 

1 U2N + U2N . ^27V+2^2Af — T2N+2T2N 



^ 4 sm^{n9) sm^{n9*) = N +-- 



2 2 2[Ti-T2] 

(45) 

Here the argument of the polynomials is 2q,/2 and for their complex conjugate it is z*/2. 
Finally, we write down finite- result for average of the absolute square of the shift, 

^ _^ 1 _ U2N + U2N ^2jv+2^2Af - T2N+2T2N 



mZa)?) = ^-^^^^ . (46) 



\D\ 

We now turn our attention to large- behavior of the result (H^ . For this purpose 
we use the ansatz fl23l) and result the flM|) for ka- Large-N behavior for the Chebyshev 
polynomials, with argument Za, may be calculated as 

exp[2iA^^„] + exp[-2iA^^„] 



T2Niza/2) 2 

n{ka] v) exp{-2ika) + [^{ka] r])r^ exp(2i/c„) ^^^^ 
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Figure 6. Asymptotic results for ( |(5iZq,P )/cr^, ( )/tT^ and 

( (3{(5iZq})^ )/ct^, shown by filled circles, squares and diamonds, vs energy index 
a. We have compared here the finite-A^ results, shown by open circles, for the exact 
resonances where N = 100 and 77 — 0.81. In the set we show these results for 14 energy 
indices near the middle of the energy band but on a different scale. 



, . , exp[ti2N + 1%] - exp[-ti2N + 1)9, 



exp{—ika) — exp{ika) 



Finally, 



T,{~z:/2) - U~zj2) = - ^l^p^ + 0{N- 



(49) 



4i 

~_cos(fcJ$5{$J, (50) 

where we have used the ansatz (l23l) in the second order polynomial ^2(2:) = 2z'^ — 1 and 
Q{<l>^} = - sin(fc^) ln{\Q\), as obtained from Eqs. ^ 

We can now plug in these results in Eq. fH6|) . These asymptotic results gives the 
numerator as [N + al + a2/{a3/N)) where al, a2/a3 are 0{N^). Similarly we obtain 
denominator as ( + N + b2) where 61 and 62 are 0{N^); see Appendix A for details. 
Thus in the leading order, we obtain 

(|(AZ,)| ) = — 1 + - — — . (51) 



What follows next is the calculation of large- A^ results for ( (JR:{AZa}y ) and 
{^{AZa})^). Since we need first to calculate {{AZa)^) and ([(AZ„)*]2), from Eq. 
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Figure 7. Shown on the same pattern of Fig. [6] but for the parametric resonances 
where N = 500 and r] = 0.81. In this figure {\5iZa\^ ) /a^ , {{^{SiZ^})"^ ) /a"^ and 
( (S{(5iZq})^ )/cr^ are shown respectively by pluses, crosses and stars. The inset is 
shown for the indices near the middle of the band. 



0]) we obtain after averaging 



and 



,El i4sin^M:) 



([(A2.)t) = ^ 



For summation we use the formula [15 



N 



n=l 



sm^fn ( 



3iV-^^(4cos[(iV+l)^^] 



TV 



sin^(n6'o 



sin(^) 

cos[2(A^ + l)g] cos(iVg) 
cos(^) 

3N - AUr,-iTN+i + 



n=l 



T2N+2U2N-I 



(52) 



(53) 



(54) 



where in the second equahty the polynomials have argument 5o/2 with 2cos(6'a) = Za- 
Similarly for the summation in Eq. ( I53l) one gets the polynomials with argument 5*/2. 
One can now use the equality fl54l) in Eqs. f l53| |53|) in order to derive finite- result 
for ( {^{AZa}y ) and ( {^{AZa}f). For large- iV we make use of the ansatz ([23]) and 
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calculate the leading order contribution as 

((a{A4}n^-|;{i-i^^^yp+.w}. (56) 

where 

~ 8 V4^ sin(A;„) {2N cos(fc„) + i sin(A;«)ln(^])) 

Equations fl5T| [55| [56]) are our main analytical results and they are also valid 
for parametric resonances with the Vt given in Eq. ( !27|) . In Fig. [6] we verify the 
asymptotic results ( ISTl l55l l56|) against their finite- counterparts, for exact resonances 
with = 100 and rj = 0.81. Fig. [7] is repeated on the same pattern but for parametric 
resonances where N = 500 and rj = 0.81. They confirm that the asymptotic results give 
a good account for the finite- results. However, there are some exception near the 
edges (not visible on the scale of the plot) where the ansatz ( l23l) is not valid. 

It turns out that in order to calculate the DOR we need the second order corrections 
{62Za), derived in Eq. ( 137|) . We have followed the method used earher [18] for Hatano- 
Nelson Model [19]. However, we have not been able to obtain a closed expression of the 
DOR. This is discussed in Appendix B where we leave the calculations with a formal 
expression for the DOR. 



j + (c.c). (57) 



VI. Numerical Methods and Verification of The Eqs. (JSH [55], [56D 

Numerical simulations for parametric resonance are always cost efficient. The reason 
being that there one deals with standard eigenvalue problem for which many fast 
subroutine packages are available, for instance LAPACK. On the other hand to verify the 
results (!5T| [55] |56|) for exact resonances, where one needs to obtain numerical solutions 
of a characteristic polynomial equation of order in a complex plane, there is no as 
good algorithm. In this paper we show results for the exact resonances by calculating 
roots of the characteristic polynomial where we have used a cost efficient numerical 
subroutine ezero. The subroutine is available on the CPC program library. There 
is one major advantage of using this subroutine over other methods, for instance the 
Newton's method. This subroutine does not require initial guesses for the roots but 
only the contour which encloses all the roots of the polynomial. Besides, it also avoids 
calculating the derivatives which may result into numerical overflow. 

In alternative to ezero we have used a different approach for calculating the 
roots. We survey the complex fc-plane for the zeros of the Det[M{k)M{ky] where 
Mrs = —2 cos{k)6rs — 'Hrs foT r,s = 1, A^ [H]. (In our system — tt < ^{k} < vr and 
'^{k} < 0.) These zeros give the eigenvalues of Ti. However, in the latter approach it 
is advisable to disintegrate the complex plane into small cells at first and then at every 
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Figure 8. Comparison of asymptotic results with numerics, for exact resonances 
where N = 100, W = 0.015 and ?/ — 0.81. In this figure, filled circles, squares and 
diamonds are the numerical results respectively for ( )/cr^, ( (3?{(5iZq})^ )/(t^ 
and ( (3{(5i-Zq})^ )/cr^ while open circles are the rescaled theories (f5T| ISS]) and (|56|) . 
In the inset we show a comparison for 14 indices near the middle of the energy band 
on a different scale of the plot. 



iteration into smaller one - only for cells which contain minima of the lowest eigenvalue 
and throwing the rest out. In this way one makes the algorithm faster and obtain the 
zeros in a reasonable precision. For a tridiagonal matrix this algorithm consumes a time 
which roughly grows with N^. However, while comparing the two methods on a simple 
machine we find that the method used in ezero is much faster than the method described 
here. We refer to [16] for further details of this subroutine. 

In Fig. IHl we compare asymptotic results with simulation done for the total number 
of realizations L = 2500, for exact resonances. In Fig. [9] we compare numerical results 
obtained for parametric resonances, where = 500, rj = 0.81 and L = 5000, with our 
theory for large- A^. Though we have considered only the flat disorder yet our results are 
valid for the Gaussian or other symmetric distribution functions. These figures show 
that our asymptotic results are in fair agreement with the numerical results for almost 
all a. For instance, near the middle of the energy band it describes reasonably well a dip 
and a peak, respectively in the ((3{AZq})^) and {{^{AZa})^). These two opposite 
effects, however, cancel out in { {\{AZa}\)^ ) ■ 

VII. Conclusion 

In conclusion, we have studied resonances in a one dimensional discrete tight-binding 
open chain in a weak disorder limit. In this study we have calculated complex energies 
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Figure 9. Shown on the same pattern of Fig. |S]but for the parametric resonances 
where N = 500, W — 0.015 and rj = 0.81. These numerical results are obtained 
from the diagonalization of TV-dimensional matrices for 5000 realizations. In this 
figure ( |(5iZq,P )/(T^, ( (5?{(5iZq})^ )/ct^ and ( (9{(5iZq})^ )/ct^ are shown respectively 
by pluses, crosses and stars while open circles are the rescaled theories ([STl [55t and 

in an open-clean chain of finite length. The result we obtain is a polynomial equation 
which we have been able to solve for long chains using an ansatz for the solution. To 
the best of our knowledge, this result has never been derived before. We have used 
a perturbation theory up to the second order where we have derived the first and the 
second order corrections to the complex energies in terms of Chebyshev polynomials. 
The first order corrections have been useful to obtain closed form of the statistical 
results for the scattered complex energies. These results have been further simplified for 
long chains where we obtain compact results. The asymptotic results have been verified 
against numerics. Evidently, in the weak disorder limit the perturbation theory predicts 
nice statistical results. Our results are new in these studies and they could be useful in 
the further studies of such systems. 

It would be interesting to study statistics of resonances in the weak disorder limit 
for higher dimensional models as well as for the cases when the site energies are not 
independent random variables but they are correlated with each other [20]. Besides, 
there has been growing interest for the case when M sites are connected to the outer 
world where 1 < M < N [2T]. We believe that our methods could be useful for the 
study of such models. Finally, we mention the case where ^(0) ~ 0{N). It requires a 
separate investigation as our perturbative analysis fails in this limit. 
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Appendix A. Large-A behavior of the denominator in (I40p 



The denominator in Eq. ( HOI) can be simplified as follows: 

D{z^) = - N + UN^i{zal2)TN+i{z^/2) - ir]e^v[I{za)] 
X ^m[k{z^)][UN-i{z^/2)f 



-N + 



-exp(-ifc^)[fi(l - r/) + 1] 



+ exp(^fc„)[^]"^(l +?7exp(2iA;„)) + 1-2?]] 



X 



-1 



Using now 
and 



2[exp(— i/ca) — exp(iA;Q,)] 

fi(l - r]) + 1 = 2 - r]exp{2ika), 

2 — 3?7 + r/^ exp(2zA;o 



n~^{l + riexp{2ika)) + 1-2// 



1 — rj exp{2ika] 



in flA.ip we get 

-N 

For I Dp this yields 



+ 



1 — rj exp{2ika) 

1 — rj cos(2A;o 



1 + r]"^ — 2ri cos(2A;q,) 



(Al) 
(A2) 



(A3) 



(A.4) 



1 + 77^ — 27/ cos(2A;q,) 

Similarly for the parametric resonances, the denominator for large N is given by 

D = - N + UN^l{Za./2)TN+l{Za/2) 

exp{—ika)^ — exp{ika)^~^ 



(A + 1/2) 
(A +1/2) 



exp{ika) — exp{—ika) 
(l + r/2)n 



2[riexp{ika) + i]"^ 



(A.5) 
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Thus for \D\ we get 



= + N 



[1 + rj'^y — Arf sm^(A;Q,) 



(A.6) 



4[1 + 7/2)2 _ 4^2 sin2(A;^)] 1 + _ sin(A;Q,) ' 
Clearly in both cases \D\'^ has a form A^^ + 51 |jvi| _|_ 52 {iVO}. 

Appendix B. Calculation of DOR 

We define the average DOR as 

where Xa = ^{Za} and = Q{Za}- Next we define Xa = ^{za}, jja = Q'l^a}, 

(SiX^) = 3fJ{(5lZ„)}, (5iy„) = 3{(5iZ„)}, (525.) = mS2Za)} and (^s^a) = 53{(52^a)} 

and then use the expansion Xa = S:a + (^i^a) + {^2X0) and 3^a = + (^i^a) + (52^0) in 
fiRT]) . We find 

TV 

{p{x,y)) = ^S{x-Xa)S{y-ya) 



N 

+ ^ [(((5i5a)((5iya))(5'(a; - Xa)5'{y - ya) - {i52Xa))S'{x - Xa)5{y - ya) 

- {{S2ya))S{x - Xa)5'{y - ya) + ]^{{hxaf)5"{x - Xa)5{y - ija) 

+ l{i6,yaf)6ix - xa)6"iy - y^)] ■ (B.2) 

Here 6'{x) = d6{x)/dx and similarly 6"{x) is the second derivative of the Dirac-Delta 
function with respect to the argument. We have already shown that ((^iXo)^) and 
{{SWa?) are of 0{a^N-^) while (((5iX«)(5iy^)) is also 0{a^N-^) since 

{{diXa){diya)) = — . (B.3) 

Motivated from pjj, we calculate the coefficient of in (1571) to obtain {{62X0)) and 
((^2^o))- Let's write 

N 

SlZa = '^enCn;a, (B.4) 

n=l 

and 

N 
n,m=l 
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Expressing the polynomials in fl4Up in terms of sinusoidal functions, we simply read-off 



Da 

cos(2n'^„) - 1 



(B.6) 



Da 

Here n' = N + 1 — n and is the denominator inside the brackets of (1401) . The 
denominator is 0{N) thus Cn-a ~ 0{N~'^), hence will be dropped off. Using now our 
ansatz, for dnn-a we obtain 

d = (1 - cos(2rv^.)) (4 .^rv^ ^ 
AN sm(Wo/A') sm(WQ) I 

- cos[(l + 2r)^J -3cos[(l -2r)^j|, (B.7) 

where r = n'/N and \1/q, = Nk^ with k^ given in f l24p . There is no further simplification 
of this result to obtain a compact and simple expression for real and imaginary parts of 
dnn;a (as the authors \18l have been able to do for with the result they obtain for the 
Hatano-Nelson model [19]). So we leave the density formally as 

N 

(p(^, y)) = poix, y)-o-^Yl {^'^^ ~ ^")^(y - y») 



a=l 

N N 

X 

n=l n=l 

where Pq{x, y) is first term of Eq. flB.2|) . 



^ ^{dnn;a} + 5{x - Xa)S'{y - y^) ^ ^{dnn;a}) , (B.8 
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